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Abstract 

A  method  for  the  modelling  of  battery  performance  under  nonstationary  load  is  presented.  For  the  purpose  of  description  of  battery  dynamic 
behaviour  a  new  equation  with  fractal  exponent  is  used.  The  applicability  of  the  proposed  battery  model  consisting  of  five  parameters,  all 
dependent  on  the  battery  state  of  charge,  is  demonstrated. 
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1.  Introduction 

The  knowledge  of  nonstationary  load  performance  is  of 
great  interest  in  the  assessment  of  the  battery  quality.  The 
difficulties  of  the  direct  experimental  measurement  of  the 
nonstationary  battery  characteristics  are  related  with  the  com¬ 
plexity  of  the  processes  taking  place  during  the  discharge. 

Most  of  the  electrochemical  methods  used  for  battery 
investigations,  such  as  measurement  of  the  volt-ampere 
(load)  characteristics,  capacity  determination,  impedance 
spectroscopy,  etc.,  gives  information  on  their  properties 
under  stationary  conditions.  Since  the  battery  state-of-charge 
is  changed  during  the  measurement  it  is  better  to  use  the  term 
quasi-stationary  instead  of  stationary  conditions. 

In  order  to  evaluate  the  evolutive  quasi-stationary  and 
dynamic  properties  of  the  battery  studied  a  special  microcycle 
test  is  applied  repeatedly  during  the  discharge  cycle  (Fig.  1 ) . 
The  microcycle  contains  two  current  steps  and  a  pause  fol¬ 
lowed  by  a  linear  sweep  segment  [  1  ] . 

2.  Experimental 

The  measurements  are  carried  out  on  a  140  Ah,  12  V  lead/ 
acid  traction  tubular-type  battery  using  computerized  test 
facility  made  in  the  authors’  laboratory,  CLEPS.  The  dis¬ 
charge  current  is  35  A.  The  test  microcycle  was  applied  at 
every  31  min  corresponding  to  13%  change  of  battery  degree 
of  discharge.  The  test  microcycle  profile  ensures  the  same 
quantity  of  electricity  as  in  the  case  of  a  constant-current 
discharge  for  the  same  time. 
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Fig.  1 .  Profile  of  the  test  microcycle. 


3.  Quasi-stationary  analysis 

The  voltage  responses  corresponding  to  the  linear  current 
sweeps  represent  a  family  of  volt-ampere  load  characteris¬ 
tics,  measured  at  continuously  progressing  state  of  discharges 
(Fig.  2) .  This  data  set  can  be  processed  by  the  technique  of 
three-dimensional  analysis  [  1  ] ,  which  takes  into  account  the 
change  of  the  state  of  discharge  during  the  measurement.  As 
a  result  the  deviation  of  s  during  the  measurement  is  avoided 
and  a  new  set  of  instantaneous  volt-ampere  characteristics  is 
produced.  These  characteristics  can  be  approximated  by  lin¬ 
ear  equations  of  the  following  form 

U(I,s)  =  U0(s)-lR(s)  (1) 

This  linear  kernel  of  the  approximation  contains  a  reduced 
open-circuit  voltage  U0(s)  and  an  effective  internal  resis¬ 
tance  /?($)>  both  function  of  s  (Fig.  3).  As  it  can  be  seen 
from  Fig.  3,  the  open-circuit  voltage,  f/0,  changes  slightly 
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Fig.  2.  Measured  volt-ampere  characteristics  at  different  states  of  discharge. 
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Fig.  3.  Dependence  of  U0  and  R  from  Eq.  ( 1 )  on  the  state  of  discharge. 

with  the  degree  of  discharge.  Its  value  decreases  with  only 
about  10%,  almost  linearly  with  s ,  while  the  effective  internal 
resistance,  /?,  depends  on  the  exponential  degree  of  discharge 
and  increases  about  two  times.  Nevertheless,  this  parameter 
could  not  be  used  directly  for  the  estimation  of  the  battery 
state-of-charge,  since  there  is  an  initial  period  of  ‘activation’ 
during  which  R  decreases.  The  simple  model,  see  Eq.  ( 1 ) ,  is 
an  effective  tool  for  modelling  and  simulation  of  the  quasi¬ 
stationary  battery  performance  and  its  evolution  with  pro¬ 
gressing  state  of  discharge. 

4.  Dynamic  analysis 
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Fig.  4.  Voltage  relaxation  after  current  step  at  different  states  of  discharge 
approximated  with  function  Fx  except  microcycle  7  fitted  with  F2. 
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Fig.  5.  Dependence  of  Rx  on  the  degree  of  discharge. 

A U(s)  =I[R{(s)  +  tf2(j)]{  1  -exp[  -  (t/T(s)n^)] }  (3) 

where  R}  (s)  represents  the  initial  voltage  step  and  R2(s)  the 
total  effective  relaxation  process. 

The  quality  of  the  fit  of  Eq.  (2)  to  the  experimental  data, 
also  processed  by  a  three-dimensional  analysis,  is  quite  good 
and  shows  stable  values  of  n  =  0.7,  changing  only  at  the  end 
of  discharge  to  n  =  0.5.  It  is  known  from  the  theory  of  topo- 
chemical  reactions  that  such  kind  of  equations  are  typical  for 
processes  limited  by  diffusion  in  the  space  with  a  variable 
dimension  [2], 

The  resistances  Rx  and  R2  also  depend  on  the  state  of 
discharge.  These  dependencies  could  be  modelled  by  empir¬ 
ical  equations  with  the  following  structure 


The  dynamic  properties  depend  also  on  the  state  of  dis¬ 
charge.  Fig.  4  represents  the  voltage  responses  of  the  current 
steps  in  the  microcycles.  In  order  to  overcome  the  difficulties 
related  to  the  classical  analysis  of  these  data  a  model  with 
fractal  exponent  is  selected 

AU(s)^Ul(s)  +  U2(s){\-txp[-(t/T(s))n^]}  (2) 

where  t  is  the  time  and  T  is  time  constant. 

Assuming  linearity,  the  voltage  transient  can  be  presented 
as 


R}  =  —  0.36634  +  4.711 32  exp (540466)  (4) 

R2=  —315.2  —  188.9s2  +  143.6 

exp(s)  4*  167.8  exp(  —s)  4-0.6597.s~1  (5) 

where  Rx  is  in  mil  and  s  is  in  parts  of  unity. 

These  dependencies  are  shown  in  Figs.  5  and  6.  The 
dependence  of  the  time  constant,  T,  on  the  state-of-charge  is 
described  by  an  empirical  equation  with  the  following 
structure 
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Fig.  6.  Dependence  of  R2  on  the  degree  of  discharge. 
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Fig.  7.  Dependence  of  the  time  constant,  T ,  on  the  degree  of  discharge. 

T=  189.1  +  1 08. 4r2  — 86.16  exp(  —5)  -97,35  exp(s)  (6) 
where  s  is  in  parts  of  unity. 

The  dependencies  of  R2  and  T  (Fig.  7)  on  s  are  quite 
complex  while  that  of  R}  is  relatively  simple.  Obviously,  the 
value  of  the  resistance  R\  could  be  used  for  estimation  of  the 
battery  state-of-charge. 


5.  Discussion 

The  knowledge  of  the  evolution  of  the  parameters  with  the 
state-of-charge  provides  for  a  simulation  of  the  battery  behav¬ 
iour  under  loads  with  different  profiles. 

The  simplest  way  of  battery  modelling  is  the  quasi-station- 
ary  approach.  Quasi-stationary  model  contains  the  relation 
between  battery  voltage  and  current  (and  as  a  result  gives  the 
open-circuit  voltage  by  Eq.  ( 1 ) )  under  stationary  conditions. 
Combining  it  with  the  so-called  ‘Ragone5  diagram  (present¬ 
ing  specific  energy  versus  specific  power)  it  is  possible  to 
determine  the  battery  abilities  to  supply  energy  under  differ¬ 
ent  load  profiles.  This  approximation  is  based  on  the  assump¬ 
tion  of  the  'additivity5  of  the  ‘Ragone5  diagram.  This 
approach  is  useful  for  relatively  fast  comparison  of  different 


Fig.  8.  Replacement  of  the  current  profile  by  a  sum  of  current  steps. 


Fig.  9.  Simulated  values  of  U  and  I  for  the  1  microcycle  under  plotted  load. 

batteries  without  long  and  expensive  experimental  work. 
Unfortunately,  this  assumption  is  not  really  true. 

One  possibility  to  avoid  this  problem  is  the  improvement 
of  the  battery  model  by  adding  a  new  dynamic  kernel.  Assum¬ 
ing  linearity  of  the  battery  behaviour  and  following  the  inte¬ 
gral  of  Duamel  cited  in  Refs.  [3,4] ,  each  load  profile  can  be 
presented  as  a  sum  of  current  steps 

I(t)=xik  (7) 

where  /k  are  step  functions. 

This  procedure  is  schematically  shown  in  Fig.  8.  The 
change  of  the  battery  voltage,  as  a  function  of  time,  caused 
by  the  kth  current  step  can  be  calculated  using  the  equation 
with  structure  similar  to  Eq.  (3) 

At/*(fej)  ^hlRi(s)  +R2{s)]  {1  -exp[  -  (tk/T(s)y^] }  (8) 

Then  the  overall  battery  voltage  is  given  by 

tf(w)=2At/*(f^)  +  t/0(5)  (9) 

where  UQ(s)  is  the  battery  open-circuit  voltage  calculated 
from  the  Eq.  ( 1 )  at  the  respective  state  of  discharge. 

Fig.  9  represents  the  simulated  voltage  and  current  behav¬ 
iour  of  the  battery  studied  discharged  by  a  power  profile 
showed  in  it  at  13%  degree  of  discharge. 

The  evolution  of  the  battery  voltage  as  a  function  of  the 
degree  of  discharge  under  the  above  load  is  shown  in  Fig.  10. 
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It  is  clear  that  the  open-circuit  voltage  ( the  first  points  of  the 
different  curves)  changed  quite  slower  than  battery  voltage 
under  some  load  (for  example  the  second  points  of  these 
curves). 

In  deep-discharge  (91%)  the  battery  still  can  follow  the 
power  profile  but  its  maximum  cause  the  battery  voltage  to 
fall  to  about  9  V.  Such  small  levels  can  cause  fast  battery 
degradation.  For  this  reason  battery  discharge  is  normally 
limited  to  about  10  V  for  six  cells  batteries.  With  this  limi¬ 
tation  the  battery  cannot  follow  the  power  profile  at  high 
degree  of  discharge  (Fig.  10) . 


6.  Conclusions 

The  battery  behaviour  under  complex  load  profiles 
depends  on  five  parameters,  each  of  them  being  a  function  of 
the  battery  degree  of  discharge 

BB[U0(s)Jtl(s)9R2(s)J(s)Ms)]  (10) 

The  above  model  takes  into  account  not  only  the  stationary 
characteristics  of  the  battery  but  its  dynamic  properties  as 
well. 

A  special  computer  software  based  on  Eq.  (9)  is  developed 
for  simulation  of  electric  vehicles  batteries  studied  under 
different  test  profiles  including  standard  ‘Euro’93’  test 
profile. 
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